UCY-PHY-94/11 



Fractal Behaviour in the 0(3) Model 

C. Christou, F.K. Diakonos, H. Panagopoulos 
Department of Natural Sciences, University of Cyprus, 
P.O. Box 537, CY-1678 Nicosia, Cyprus 



Abstract 

We study domain formation in the two-dimensional 0(3) model near criticality. The 
fractal dimension of these domains is determined with good statistical accuracy. 



1 Introduction 



The D = 2 nonlinear 0(3) sigma model is known to be a prototype model for asymptotically 
free theories like the D = 4 non-abelian gauge theories. In previous studies of SU{2) lattice 
gauge theory [1] it has been shown that the domains of deconfinement (clusters with vanishing 
values of the Polyakov loop) have, at the critical temperature, a fractal structure with dimension 
Dp 2.5. Given the close similarities between the SU{2) gauge theory and the 0(3) nonlinear 
sigma model, one can naturally ask whether such fractal structures can also emerge in the 0(3) 
model. There are however some features of the 0(3) lattice theory which make the search for 
such fractal domains a difficult task : 

1) There is no phase transition at finite temperature. One can of course consider T = 
as the critical temperature because the correlation length diverges as T tends to zero. On the 
other hand this means that the order parameter will vanish at any finite temperature. A similar 
situation occurs in the 1 — d Ising model. In that case, a study of intermittency indices [2], 
possibly connected with an underlying fractal geometry, shows that even in the presence of an 
external field, where the phase transition occurs at finite temperature, there are no self-similar 
structures with fractal geometry surviving in the thermodynamic limit. 

2) As T tends to zero, one observes large domains each having a uniform orientation of the 
spins and the system near the critical point tends to magnetize spontaneously. The underly- 
ing rotational symmetry of the model is expressed however through long-wavelength fluctua- 
tions corresponding to a slowly rotating magnetization as well as short-wavelength fluctuations 
around this slow evolution. It is the presence of these fluctuations which hinders the definition 
of the ordered domains. 

Using the so-called Symanzik tree improved action we have studied on the lattice the for- 
mation of domains with spontaneous magnetization near the critical temperature. We consider 
here the situation where the lattice size L is of the order of few correlation lengths. A suitable 
order parameter is introduced and the dimension of domains with non- vanishing value of the 
order parameter is determined. It turns out that the geometry of these domains is described 
by a fractal dimension Dp 1.75 which is configuration independent. 

The paper is organized as follows : In section 2 we give the action of the model considered 
and the definitions of the observables used in the subsequent calculations. In section 3 we de- 
scribe the algorithm used for the determination of the observables and of the fractal dimension. 
Finally in section 4 we present our numerical results and we discuss their consequences. 

2 Lattice formulation of the model 

We use the lattice action : 

(3 4 1 

S sym = "^El^nC+f. ~ (!) 

where cf> a is a 3-component scalar field subject to the nonlinear constraint : <f> a (p a = 1 at each 
lattice site. The model is a special combination of a nearest and a next-to-nearest neighbour 
action chosen so that the short distance (ultraviolet) behaviour is markedly more similar to the 
behaviour of the continuum action [3]. The correlation length £ defined through the mass gap 
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ma [4] : £ = I /ma diverges exponentially as the temperature tends to zero. Near the critical 
point, configurations with domains of non vanishing magnetization occur with radii varying 
from the lattice spacing up to the magnitude of the correlation length. The mean magnetization 
in each independent configuration has a different orientation. Due to the rotational symmetry 
of the action configurations with different orientation of the mean magnetization are equally 
probable. We characterize each configuration by the magnitude M of its mean magnetization: 

\r jzY.° : - w (2) 

L n 

After equilibrium, M is statistically distinct from that of random configurations; the latter 
typically lead to values of M smaller by 2 orders of magnitude (see Fig.f). This suggests the 
introduction of the constraint effective potential (CEP) U e ff [5] which is a function of M only 
and is defined by rewriting the finite lattice path integral in the form : 

Z = j ,l.\l .\l-< 1 •' :Wi (3) 

The free energy differences for various configurations of the theory are proportional to the 
difference of the corresponding values of the CEP only. The form of CEP is similar to that of 
the effective potential of a theory undergoing a second-order phase transition. We conjecture 
that this property is the indication of the self-similarity in the geometry of the ordered domains 
which in turn leads to the fractality of the system. 

Given a configuration with spontaneous mean magnetization of magnitude M , its domains 
are determined by evaluating the quantity : 

(b a M a 

P = 0(^_ _ e ) ( 4 ) 

P plays the role of the order parameter in our considerations. We introduce also a cut-off 
parameter e in order to take into account the precision of the spins due to small wavelength 
fluctuations. Varying e we can estimate the effect of these fluctuations to the calculated fractal 
dimension. If an underlying fractal geometry is present we expect Dp to depend very weakly 
on e. 



3 Determining the fractal dimension 

After equilibrium is achieved we use Eq.(4) to determine the ordered clusters for each indepen- 
dent configuration. The resulting geometry of the clusters for some configurations is shown in 
Fig. 2. For each configuration we calculate the fractal dimension of the ordered domains. We 
define the diameter of each domain D through 

dp = max |x; — xj|, Vx;, Xj £ D (5) 

and we include in our calculation only domains greater than o? m i n , where o? m i n is determined as 
a fraction of the correlation length £, such that a) the results are not affected by the presence 
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of pointlike domains and b) we consider clusters of many different sizes. In actual runs we took 
^min = £/15- The fractal dimension is determined using the correlation function defined in [6] : 

I N 

C ( R ) = ^iE 0(i2-|xi-xj|) (6) 

* J=i 

where the summation extends over all points with P = 1 (cf. Eq.(4)). For a wide range of 
values of R one expects: 

C(R) = aR Dp (7) 

Indeed, our data exhibit this behaviour, for d m i n <R < L/2. A typical sample is shown in Fig. 
3. Fitting our data to the power law (7) we extract Dp with a statistical error coming from an 
average over configurations. It is clear from the figure that changing e has a negligible effect 
on Dp. The calculated fractal dimension shows a weak dependence on the configuration. 

It is interesting to study the dependence of the fractal dimension on the temperature. By 
performing the calculation of Dp at a higher f3 we must take into account the scaling properties 
of the system. Therefore, due to the increase of £, which is the only existing physical length 
scale, we must renormalize the value of e so that the domains which we consider for determining 
Dp remain of the same size in units of £. One way to do this, is to calculate the mean diameter 
of the domains for all configurations, defined by : 
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where the summation is over the domains D and the configurations C . Np> is the number of 
points in each such domain. We then determine e by requiring: 

d . d. 

Notice that the size of the domains is a function of e. We find that the geometrical picture 
produced with the adjusted value of e leads us to approximately the same value of Dp. 



4 Numerical results and conclusions 

Our production runs were carried out at 3 different values of /3, listed in Table I. We used 
the "overheat bath" updating algorithm [7] to thermalize our configurations. The exponential 
autocorrelation time r exp [8] was found to be < 1000 for all /3's used. After a few thousand ther- 
malization sweeps we used a mixture of microcanonical(80%)-overheat bath(20%) algorithms 
for our measurements. Table I also presents the autocorrelation times T 8nt , using the mixed 
algorithm, for the various /3-values. For each f3 we produced an ensemble of 50 independent 
thermalized configurations, measuring every 50 sweeps for the first two f3- values and every 100 
sweeps for f3 = 1.68, which was sufficient to determine Dp with a statistical error smaller than 
1%. 

For any finite value of /3, one may expect to observe fractal behaviour at best up to some 
maximum length scale: So long as £ is finite, it is easy to see that, for R ^> £, definition (6) 
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will lead to C(R) oc R De , where D e is the embedding dimension (D e = 2 in our case). Thus, 
further data (at larger R) in Fig. 3 will necessarily tend toward a line with slope 2. To find the 
range of length scales with fractal behaviour, one must determine the greatest interval of values 
for i?, taken around £, such that the corresponding data in Fig. 3 can be fitted by a straight 
line at a certain confidence level. Clearly, lattice sizes L greater than this maximum range 
are unnecessary. We chose L so that L/£ is a constant, for uniformity in comparing results at 
different f3. We took L/£ = 5; the corresponding values are listed in Table I. With this choice 
of L, fractality was seen to extend out to at least L/2. The values of £ given in Table I are 
taken from [4]. Any error in the measurement of £ can influence the precision in the estimation 
of d/£, but has little effect on the measured value of the fractal dimension. To establish a 
fractal dimension in this system it is imperative to verify its independence from e, in view of 
the fact that domain size is e-dependent. The allowed range of values for e is dictated by the 
average domain size, which must be neither pointlike nor comparable to the lattice itself. We 
verified that, within this range, Dp is indeed e-independent. We list in Table I the results for 
two sets of values for e, chosen in a way as to satisfy Eq.(9) at different f3. The errors reported 
in Table I are purely statistical, while a comparison at different f3 gives an estimate of possible 
systematic errors. 

The study presented in this paper can be taken over to simulations of QCD with different 
fermionic content. These afford us with a rich spectrum of phase transitions, some of which 
could show fractal behaviour. It would also be interesting to compare the fractal dimension in 
these cases to the critical indices coming from intermittency, as observed in heavy ion collisions. 
We expect to return to this subject in the future. 
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Figure captions 



Fig. 1 : Distribution of M over different configurations at f3 = 1.68. The solid spike is the 
distribution corresponding to random configurations. 

Fig. 2 : A typical pattern of domains. (/3 = 1.68, e = 0.8) 

Fig. 3 : C(R) vs. R for different values of e. 



Table 1 





Tint 

(sweeps) 




L 


e 


dH 


D F 


1.44 


24 


19.9 ± 0.2 


99 


0.850 
0.900 


0.198 ± 0.006 
0.142 ± 0.003 


1.742 ± 0.011 
1.740 ± 0.018 


1.56 


38 


41.7 ± 0.3 


210 


0.790 
0.845 


0.195 ± 0.004 
0.138 ± 0.002 


1.760 ± 0.012 
1.753 ± 0.012 


1.68 


75 


83 ± 2 


399 


0.750 
0.800 


0.195 ± 0.006 
0.139 ± 0.004 


1.779 ± 0.010 
1.785 ± 0.014 
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